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We present the creation and time evolution of two-dimensional Skyrmion excitations in an antifer- 
romagnetic spinor Bose-Einstein condensate. Using a spin rotation method, the Skyrmion spin tex- 
tures were imprinted on a sodium condensate in a polar phase, where the two-dimensional Skyrmion 
is topologically protected. The Skyrmion was observed to be stable on a short time scale of a few 
tens of ms but to have dynamical instability to deform its shape and eventually decay to a uni- 
form spin texture. The deformed spin textures reveal that the decay dynamics involves breaking 
the polar phase inside the condensate without having topological charge density flow through the 
boundary of the finite-sized sample. We discuss the possible formation of half-quantum vortices in 
the deformation process. 
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Skyrmions are particle-like topological solitons, first 
proposed to account for the existence of protons and 
neutrons in high energy physics [lj and recognized to 
have an important role in many condensed matter sys- 
tems such as liquid 3 He-A [2 , quantum Hall systems [3], 
liquid crystals [4], and helical ferromagnets [5 , 6 . Atomic 
Bose-Einstein condensates (BECs) with internal spin de- 
grees of freedom are attractive for the study of topologi- 
cal objects. The rich structure of their order parameters 
can accommodate various topological excitations [7HT2]. 
Furthermore, precise spin manipulation techniques have 
been developed to prepare topological spin structures of 
interest, providing unique opportunities to study their 
stability and dynamics [T3HT9] . 

Stability is one of the key issues in the study of 
Skyrmions [8] I20H22] . 3D Skyrmions, which are not ex- 
perimentally discovered yet, are anticipated to be en- 
ergetically unstable against shrinking to zero size with- 
out additional stabilizing mechanisms. In the search for 
stable Skyrmions, spinor BECs are particularly appeal- 
ing because of the possibility of engineering experimen- 
tal conditions such as interaction properties and exter- 
nal potentials including rotation [21] . Earlier experimen- 
tal efforts have been focused on creating 2D Skyrmion 
spin textures in spinor BECs by using phase imprinting 
methods [15] [19], where, however, the BECs were not 
trapped so that the stability of the spin textures could 
not be investigated. It is important to note that the fer- 
romagnetic phases with SO (3) symmetry studied in the 
previous experiments [I5j [19] have the trivial second ho- 
motopy group, meaning that the spin textures are not 
topologically protected [23] . The coreless vortex states 
with 2D Skyrmion spin textures, referred as Anderson- 
Toulouse [24] or Mermin-Ho vortices [25 , have infinite 
energy in a 2D plane and therefore fundamentally dif- 
ferent from the Skyrmions, localized topological solitons 
with finite energy [2T] . 



In this Letter, we demonstrate the creation of 2D 
Skyrmions in an antiferromagnetic F = 1 spinor BEC 
in a polar phase and study their time evolution in a 
harmonic potential. The order parameter manifold of 
the polar phase is M = (U(l) x S ,2 )/Z 2 , having the 
nontrivial second homotopy group ^(M) = Z [23| [26]. 
Thus, the spin textures realized in this work represent the 
first example of topologically stable 2D Skyrmions in a 
spinor BEC. Their nature is confirmed by the observation 
that the spin textures are stable on a short time scale. 
They are, however, dynamically unstable to deform and 
eventually decay to a uniform texture. The spatial pat- 
terns of the deformed spin textures reveal that the decay 
dynamics involves breaking the polar phase inside the 
Skyrmion without having topological charge density flow 
via the boundary of the condensate. We suggest that 
half-quantum vortices might be nucleated with the non- 
polar defects in the deformation process. 

In terms of a unit spin vector <i, the 2D Skyrmion spin 
texture with z-axis symmetry is given as 

d(r, (j)) = cos f3(r)z + sin j3{r)r (1) 

with the boundary conditions, /3(0) = and /3(oo) = tt. 
This spin texture has the topological charge 

Q = J dxdy d- (d x dx d y d) = 1 (2) 

which represents the number of times the spin texture 
encloses the whole spin space. The order parameter of 
a F = 1 BEC can be written as \I> = (^i, ^-i) T = 
y/ne L '®( > , where ip m is the \m z = m) component of the or- 
der parameter (m = 0, ±1), n is the atomic number den- 
sity, $ is the superfluid phase, and ( is a three-component 
spinor. The ground state of an antiferromagnetic BEC 
is polar, i.e. |rap = 0) [27]. Setting d as the spin 
quantization axis, ( = ( ~ dx + ld y ? ^ d x +id y anc j ^ e 2D 
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FIG. 1: (color online) Illustration of the creation process of 
a 2D Skyrmion in a spinor Bose-Einstein condensate, (a) A 
quasi-2D condensate is prepared in an optical trap and a 3D 
quadrupole magnetic field is generated from a pair of coils 
with opposite currents. An additional pair of coils provides 
a uniform bias field to control the z-position of the zero-field 
center of the quadrupole field, (b) The zero- field point pene- 
trates through the condensate and the atomic spin is driven to 
rotate by the change of the local magnetic field direction. The 
tilt angle of the atomic spin depends on its radial position, 
resulting in a 2D Skyrmion spin texture. 



Skyrmion spin texture in a polar BEC is 
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The spin winding structure around the z-axis is expressed 
with the opposite phase windings of the \m z = 1) and 
\m z = —1) components. Note that the total angular 
momentum of this spin texture is zero. 

For creating a 2D Skyrmion spin texture, we employed 
a spin rotation method [15][T6J[28] as illustrated in Fig [I] 
A quasi-2D, polar BEC confined in an optical dipole trap 
is placed on the z=0 plane and a 3D quadrupole magnetic 
field is generated as 



B(r, z) = B'(rr - 2zz) + B z z, 



(4) 



where B' is the radial magnetic field gradient. The z- 
position of the zero-field center of the quadrupole field is 
controlled by the axial bias field B z . Initially, the spin 
vector of the polar BEC is d= +z with B z > B' R > 0, 
where R is the radial extent of the conensate. By increas- 
ing B z to B z <C —\B'R\ < 0, the zero-field point pen- 
etrates through the condensate and the magnetic field 
on the condensate rotates by tt. In the outer region 
R > r > r c ~ [h\B z \/ii B B ,2 ) 1/2 , where ji B is the Bohr 
magneton, and H is the Planck constant divided by 27r, 
the local field direction changes slowly with respect to 
the local Lamor frequency and d(r) adiabatically follows 
the local field direction to — z. On the other hand, in the 
center region r <C r c , the local field rotates so abruptly 
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FIG. 2: (color online) Creation of 2D Skyrmions in spin-1 po- 
lar Bose-Einstein condensates, (a) — (c) Density distribution 
of the \m z = 0, =bl) spin components were measured by tak- 
ing absorption images after a Stern-Gerlach separation. The 
field ramp rate was (a) \B Z \ =12, (b) 32, and (c) 80 G/ms 
with a quadrupole field gradient B' =8.1 G/cm. The contour 
plots (d) — (f) show the distribution of the tilt angle /3(x,y) 
in (a) — (d), respectively. (g,h) The density profiles of the 
spin components are the horizontal center cuts in (a,c). The 
Skyrmion size R^/2 versus (i) \B Z \ and (j) B'' 1 . The solid 
lines are square-root and linear fits to the data points, respec- 
tively. The field of view in (a) — (c) is 1.2 mmx 330 /xm. 



that d{r) cannot follow the field direction. It is obvious 
that d(0) keeps its direction in +z. Consequently, the tilt 
angle f3{r) of d(r) continuously changes from /3(0) = 
to P{R) = 7r, satisfying the boundary conditions of the 
Skyrmion spin texture. 

Bose-Einstein condensates of 23 Na atoms were gener- 
ated in the \F — l^mp — —1) state in an optically 
plugged magnetic quadrupole trap [29] and transferred 
into an optical dipole trap formed by focusing a 1064-nm 
laser beam with a 1/e 2 beam waist of 1.9 mm(17 /im) 
in the y(z)-direction. Further evaporation cooling was 
applied by lowering the trap depth and a quasi-pure con- 
densate of 1.2 x 10 6 atoms was obtained. The conden- 
sate was prepared in the \m z = 0) state by using an 
adiabatic Landau- Zener RF sweep at a uniform bias field 
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FIG. 3: Interference of the coreless spin vortex state, (a) 
Initial density distribution of the three spin components, (b) 
A RF pulse was applied before a Stern- Gerlach separation 
and fork-shaped interference fringes appear in the \m z = 0) 
atom cloud, (c) Numerical simulation of our experimental 
condition with the phase winding numbers ( — 1,0, 1) for the 
\m z = +1), |0), and | — 1) components. 



B z = 21 G. The trapping frequencies of the final optical 
trap were {uo x , uo yi uj z )=2ttx (3.5, 4.6, 430) Hz and the 
transverse Thomas-Fermi radii of the trapped conden- 
sate were (R x , R y )~(150, 120) jam. For a typical atom 
density n = 1.2 x 10 13 cm -3 , the spin healing length 
£ s « 40 /im [30] which is much larger than the thickness 
of the condensate, so the spin dynamics in the condensate 
is of 2D character. 

For the spin texture imprinting, the quadrupole field 
was adiabatically turned on to B' = 8.1 G/cm in 40 ms 
with B z = 500 mG and the the axial bias field was 
rapidly ramped to B z < —500 mG at a variable ramp 
rate B z . Then the quadrupole field was switched off 
within 150 /is and the axial bias field was stabilized to 
B z = —500 mG within 100 /is. The structure of the 
spin texture was determined by measuring the density 
distributions no,±i(#>2/) of the \m z = 0, ±1) components 
after a Stern-Gerlach spin separation. After switching 
off the optical trap, the magnetic field was adiabatically 
rotated and a field gradient was applied for 5 ms in the x- 
direction to spatially separate the three spin components. 
Then atoms are pumped into the \F = 2) state and an 
absorption image was taken using the \F = 2) — >> \F' = 3) 
cycling transition. During the total 15-ms time-of-flight 
the condensate expanded transversely by less than 10% 
and we assume that the measured density distributions 
adequately reveal the in-situ spin texture. 

Ring-shaped spin textures were observed in the con- 
densate after the imprinting process (Fig. [2). The |0) 
component shows a clear density-depleted ring and the 
| ± 1) components occupy the ring region with equal 
densities. The radius of the ring, characterizing 

the size of the spin texture, could be controlled with 
the field ramp rate B z and the quadrupole field gra- 
dient B' \ as R n /2 oc {Bzl 1 / 2 / B f , having the same de- 
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FIG. 4: Deformation and decaying of the 2D Skyrmion in a 
harmonic potential. The time evolution of the spin texture 
was measured with a variable dwell time in an optical trap. 
For this data, \B Z \ =20 G/ms and B' =8.1 G/cm. The field 
of view is 1.0 mmx 330 /im. 



pendence of the critical radius r c of the adiabatic spin 
rotation. This validates the aforementioned qualitative 
description of the imprinting process. The tilt angle 
/3(x,y) was reconstructed from no,±i(x,y) with the rela- 
tion cos 2 f3 = no/(ni + no + n_i). When R^/2 < 0.3ir! x , 
| cos f3(R x )\ > 0.95, satisfying the boundary condition of 
the 2D Skyrmion within our imaging resolution. When 
Rn/2 ~ Rxi the outer ring part of the |0) component dis- 
appears and a coreless spin vortex state is formed, where 
the stationary |0) atoms fill the core of the | ± 1) vortices 
with opposite circulations [32] 

The spin winding structure of the spin texture was 
confirmed from the matter wave interference between the 
spin components. A 14- /is RF pulse of 430 kHz was ap- 
plied at £^=620 mG before turning off the optical trap. 
The power and duration of the RF pulse were set with a 
\m z =0) condensate to transfer all atoms equally into 
the \m z — ±1) states. Three-to-one, fork-shaped in- 
terference fringe patterns were observed (Fig. [3|, clearly 
demonstrating that the relative phase winding between 
the | zb 1) components around the core |0) component is 
47r [33[ [34]. Since the total angular momentum of the 
condensate should be conserved to be zero in the im- 
printing process [35], the phase winding numbers of the 
spin components are ( — 1, 0, 1). Together with the radial 
distribution of /3, this clearly shows that a 2D Skyrmion 
is indeed created in the polar condensate. 

In order to study the stability of the 2D Skyrmion, 
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we measured the time evolution of the spin texture with 
a variable dwell time t in the optical trap and observed 
that it dynamically deforms and decays to a uniform spin 
texture (Fig. [4| . Modulation of the ring structure gradu- 
ally developed and distorted the whole spin texture over 
100 ms. Eventually, the spin texture became uniform in 
the \m z =0) state at t > 300 ms. Recent theoretical cal- 
culations have shown that 2D Skyrmion excitations are 
energetically unstable to decay to a uniform spin texture 
by expanding or shrinking in a harmonic potential [22] 
but we estimate that the thermal dissipation is negligible 
in our experiments for the upper bound of the conden- 
sate temperature < 65 nK. The decay of the Skyrmion is 
a dynamically driven process. 

It is remarkable that the outer |0) component kept 
surrounding the | ± 1) components in the decay dynam- 
ics, implying that there was no topological charge den- 
sity flow via the boundary of the finite-sized conden- 
sate. Since it is topologically impossible for the isolated 
Skyrmion to unwind to a uniform spin texture when the 
global polar phase is preserved, non-polar defects must 
have developed inside the condensates during the decay 
dynamics. Breaking the polar phase might be attributed 
to the quadratic Zeeman effects and the induced spin cur- 
rents [10 . The spin relaxation time of a \m z = ±1) equal 
mixture was measured to be over 1 s at B z = 500 mG, 
clearly indicating that the rapid development of non- 
polar defects is due to the presence of the Skyrmion spin 
texture. 

One interesting observation in the deformed spin tex- 
tures is that the core |0) component could be connected 
to the outer component or divided into two parts (Fig. [5]). 
Because of the symmetry of the order parameter of the 
polar BEC under (d, #) — >• (— d, i? + 7r), the two regions 
with a 7r disclination in d can be continuously connected 
via a 7r change of the superfluid phase. This observa- 
tion suggests possible formation of half-quantum vortices 
around non-polar local defects in the deformation process 
as described in Fig. [5] The | ± 1) components seemed 
to keep occupying the density-depleted region of the |0) 
component without developing noticeable local spin po- 
larization 774 — W-ij supporting the possibility. A similar 
situation has been studied with a point defect in a po- 
lar 3D BEC, where the 't Hooft-Polyakov monopole con- 
tinuously deforms to a half-quantum vortex ring (Alice 
ring) [[36] . 

In conclusion, we have demonstrated the creation of 
the topologically stable 2D Skyrmions in a polar BEC 
and presented the first measurement of the dynamics 
of the Skyrmion in a harmonic potential. This work 
opens a new perspective for the study of Skyrmion dy- 
namics in spinor BECs. The experimental technique can 
be extended to create multiple Skyrmions in a conden- 
sate or other exotic spin textures such as knots [10] and 
monopoles [36] in 3D condensates. 

The authors thank W. Ketterle and J. H. Han for dis- 
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FIG. 5: Dynamic formation of a pair of half-quantum vor- 
tices. The line segments represent the axis of the spin vector 
d and the grey area indicates the region where the \m z = 0) 
component mainly occupies. The non-polar cores of the half- 
quantum vortices are marked by and x , implying opposite 
windings of the superfluid phase. Examples of the dynamic 
evolution of the 2D Skyrmion spin texture at (c) t = 30 and 
(d) 50 ms, displaying density profiles similar to those in (a) 
and (b). (e,f) Same as (c,d) with guidelines for the peak den- 
sity regions in the \m z = 0) atom cloud. 
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